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We use diagrammatic many-body perturbation theory in combination with low-momentum inter- 
actions derived from chiral effective field theory to construct effective shell-model transition oper- 
ators for the neutrinoless double-beta decay of '^^Ge and *^Se. We include all unfolded diagrams 
to first- and second-order in the interaction and all singly folded diagrams that can be constructed 
from them. The resulting effective operator, which accounts for physics outside the shell-model 
space, increases the nuclear matrix element by about 20% in ^®Ge and 30% in *^Se. 



PACS numbers: 23.40.-s, 21.60. Cs, 24.10. Cn, 27.50.-|-e 



I. INTRODUCTION 

The experimental discovery of neutrinoless double- 
beta (Oi^/3/3) decay, a nuclear-weak process that occurs 
extremely slowly if at all, would have deep implications 
for particle physics. Since Qv^P decay can occur only 
if the neutrino is its own antiparticle, an observation 
would at once establish the neutrino as a Majorana par- 
ticle. Furthermore, from a measured lifetime we could, in 
the absence of exotic new physics, determine an average 
neutrino mass rrii, = ^^ Ueimf, where i labels the mass 
eigenstates, and U is the neutrino mixing matrix [T]. This 
mass cannot be extracted from a lifetime, however, with- 
out first knowing the value of a nuclear matrix element 
that also plays a role in the decay. The entanglement 
of nuclear and neutrino physics has led to a small but 
concentrated effort within the nuclear structure commu- 
nity to calculate the nuclear matrix elements, which are 
not themselves observable. While various theoretical ap- 
proaches agree to within factors of two or three, — a 
range many structure theorists might find not unreason- 
able — the uncertainty in the effective mass that can 
be extracted from an observed lifetime is at least that 
large as a result. Since large-scale experiments will be 
reporting results in the coming years, we need to work 
quickly to improve the accuracy of the matrix-element 
calculations. 

Of the theoretical methods currently employed, the nu- 
clear shell model is the only approach that offers an ex- 
act treatment of many-body correlations, albeit within 
a truncated single-particle (valence) space above some 
assumed inert core. Though most of the physics govern- 
ing double-beta (/3/3) decay indeed resides in this valence 
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space, correlations involving neglected single-particle or- 
bitals may contribute non-negligibly to both the Hamil- 
tonian and the Oz^/3/3-decay transition operator, each 
of which is a basic ingredient in any nuclear matrix- 
element calculation. Contributions to the Hamiltonian 
can be, and have been, included in the construction of 
an effective valence-space Hamiltonian, TJoff , through dia- 
grammatic many-body perturbation theory (MBPT) [5] . 
But the analogous contributions to an effective valence- 
space 0^'/3/3-decay operator, with the exception of a crude 
renormalization of gA, have thus far been almost com- 
pletely ignored. The first and only work to apply MBPT 
to the Oj^/3/3-decay operator considered only diagrams 
that were first order in the interaction (a G-matrix) , plus 
a few selected higher-order contributions [3]. 

In this article we carry out a much more compre- 
hensive computation, providing the first steps towards 
a true first-principles calculation of nuclear matrix ele- 
ments based on chiral nuclear forces [4l. We first define 
and compute an X-box consisting of all diagrams now 
to second order in the interaction and in a much larger 
Hilbert space than used in Ref. [3] . We then consider con- 
tributions of folded diagrams together with state norms, 
which must be explicitly computed for effective transition 
operators (see Eq. Q). We finally apply the resulting 
two-body effective Oi^/3/3-decay operator, together with 
wavefunctions from existing shell-model calculations, to 
obtain corrected nuclear matrix elements in the p/-shell 
Oi^/3/3-decay candidates ^^Ge and ^^Se. 

Assessing the accuracy of the perturbative expansion 
is a central challenge for MBPT. Although a nonpertur- 
bative treatment of core polarization found only mod- 
est changes in if off [Sj , the analogous impact on effective 
/3/3-decay operators is unclear, and ultimately a nonper- 
turbative method that goes beyond core polarization, al- 
lowing controlled approximations to both the effective 
Hamiltonian and transition operators, will be preferable. 
Coupled-cluster theory |S1 [7] and the in-medium similar- 



ity renormalization group [8HID] are promising nonper- 
turbative methods, but neither has yet been applied to 
/3/3 decay. The situation may be different in a few years, 
but at present MBPT is still the best method to inves- 
tigate microscopic many-body corrections to the shell- 
model 0:^/3/3-decay operator. And even within MBPT, 
as we have noted, there is essentially no work, outside 
of Ref . [3j, on two-body transition operators, making the 
topic almost completely unexplored. 

The remainder of this paper is structured as follows: 
Section |ll] describes the ingredients of our calculation, 
including definitions of the matrix elements we compute, 
the framework for obtaining the nuclear interactions with 
which we begin, and the details of our many-body formal- 
ism for calculating effective /3/3-decay operators. Section 
III presents our results for ^^Ge and *^Se, updating the 
matrix element for *^Se first reported in Ref. [H]. Fi- 
nally, Section |IV| discusses the significance of the results 
and outlines steps that will improve their accuracy. 



II. METHODS 
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and rup denotes the proton mass and TO,r the pion mass. 
The closure approximation is not good for two- 
neutrino double-beta {2vj3/3) decay, which we briefiy dis- 
cuss later. The matrix element governing that process 
contains a complete set of intermediate states, viz.: 



M: 



2u 



Y^ {f\J2a^aTa 



I Efc ^bT,^ 



En - (M, + Af/)/2 



(4) 



where n denotes states in the intermediate nucleus with 
energy En , Mi and Mf are the masses of the initial and 
final nuclei, and the effects we neglect (e.g., forbidden 
currents, the Fermi matrix element, etc.) are small. We 
are unable to obtain a complete set of intermediate states, 
so we can treat 2z//3/3 decay only in the closure approxi- 
mation, viz.: 
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A. Decay Operator 

In the closure approximation (which is good to at worst 
10% or so [E]), the nuclear matrix element governing 
Oi//3/? decay can be represented as the ground-state-to- 
ground-state matrix element of a two-body operator. Ne- 
glecting the so-called "tensor term," the effect of which 
is a few percent [131 [T4] , the matrix element is given by 
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where \i) and |/) are the ground states of the initial and 
final nuclei, rab = \fa — T^bl is the distance between nu- 
cleons a and b, jo is the usual spherical Bessel function, 
and the nuclear radius R is inserted to make the matrix 
element dimensionless, with a compensating factor in the 
phase-space integral that multiplies the matrix element. 
The "form factors" hp and hcT Eire given by 
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where 
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Although the approximation is poor, and we cannot use 
it to deduce the real 2z//3/3-decay matrix element, the clo- 
sure matrix element and the real one change in a similar 
way when correlations are added. 



B. Nuclear Interactions 

Diagrammatic MBPT was reviewed extensively some 
years ago [T3 [TC], but since then, driven by ad- 
vances in chiral effective field theory (EFT) [4j and 
renormalization-group (RG) methods |T7j, it has seen 
something of a revival [IH]. Chiral EFT is a systematic 
expansion of nuclear interactions and electroweak cur- 
rents in which three- (3N) and higher-body forces arise 
naturally. Beginning from the chiral two-nucleon (NN) 
potential of Ref. ^19J . we construct a low-momentum 
interaction (T^owfe), with cutoff A = 2.0 fm^ , via RG 
evolution '17' '20' , explicitly decoupling high-momentum 
components from those at the nuclear-structure scale. In 
contrast, the G-matrix [TS], often taken as a starting 
point in nuclear structure calculations, deals with high- 
momentum modes by particle-ladder resummation, and 
does not adequately decouple low- from high-momentum 
degrees of freedom. As a result, many-body methods 
based on Viowfc tend to converge better than those us- 
ing a G-matrix 6 . Recent work with MBPT based on 
low-momentum NN-I-3N interactions has led to the de- 
velopment of non-empirical valence-space Hamiltonians 
for proton- and neutron-rich systems [2TH25] . While 3N 
forces are neglected here, we plan to include them in our 
future 0j//3/3-decay nuclear-matrix-element calculations. 

The one drawback of using low-momentum interac- 
tions in calculations of effective operators is that high- 
momentum physics cannot be included explicitly. The 
effects of high-momentum (short-range) correlations on 



the 0!//3/3-decay operator are both small and now well 
understood, however, and we include them via an effec- 
tive Jastrow function that has been fit to the results of 
Brueckner-theory calculations [55]. 



C. Effective Two-Body Transition Operators 

As we have noted, existing work on MBPT contains 
little about effective two-body operators other than the 
Hamiltonian, where Refs. [2J H?! [5S] provide the most 
comprehensive discussion. No matter the two-body op- 
erator of interest, however, the starting point is always 
the construction of projection operators P and Q that di- 
vide the full many-body Hilbert space into a model space, 
in which subsequent exact diagonalization is carried out, 
and everything else. In our calculations in nuclei with 
mass near A = 80, the model space consists of the O/5/2, 
lP3/2: 1^1/2 7 and Ogg/2 single-particle orbits, for both 
protons and neutrons, above a ^^Ni core in a harmonic- 
oscillator basis of 13 major shells with fnu = 10.0 MeV. 

After specifying the model space, one must define a 
mapping between eigenstates of the full Hamiltonian and 
projections of those eigenstates onto the model space. In 
MBPT this is done perturbatively. The result is a set of 
diagrams with two incoming legs and two outgoing legs, 
with each diagram representing a contribution to the two- 
body matrix elements of the effective Hamiltonian or ef- 
fective (two-body) transition operator. The usual Feyn- 
man rules are used to evaluate the diagrams, but to the 
set of familiar-looking diagrams one must add "folded" 
diagrams, which eliminate the energy dependence of the 
effective operator [15l [16] . One way to organize the sum 
of all diagrams is by grouping all those without folds into 
a "Q-box" (for the Hamiltonian) or an "X-box" (for the 
transition operator) and then writing the complete sum, 
including folded diagrams, in terms of the Q- and X- 
boxes and their derivatives with respect to unperturbed 
energies. The first few terms in the Q- and X-boxes ap- 
pear in Figs. [1] and [2] 

Folding is significantly more complicated for a two- 
body transition operator, which combines X- and Q- 
boxes, than for the Hamiltonian, where only Q-boxes are 
needed. Effective model-space operators in the basis of 
energy eigenstates are always defined (for a bare operator 
M) via 
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where the states that lie in the model space, \ics) = P \i) 
and |/ef=f) = P\f)j are not in general normalized. If M is 
the Hamiltonian, then only diagonal matrix elements are 
nonzero, and the denominator is canceled by a similar 
factor in the numerator. For two-body transition opera- 
tors, that is not the case, and state norms must be explic- 
itly computed. Prior authors have approached the issue 
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FIG. 1. The Q-box to second order in V^ow k (ellipses indicate 
higher-order terms). The first line contains one-body con- 
tributions and the others two-body contributions. Exchange 
diagrams, though not shown, are included in our calculations. 



of norms in several ways. References [2] and [28], for in- 
stance, choose to expand the denominators and fold them 
into the numerators, thus completely eliminating all dis- 
connected diagrams. The resulting expressions, however, 
become complicated as the number of folds increases, and 
the approach requires the construction of a special basis 
as an intermediate step. For these reasons Ref. [571 advo- 
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FIG. 2. (Color online) The X-hax to first order in Viowfc. 
Solid (red) up- or down-going lines indicate neutrons and dot- 
ted (blue) lines protons. The wavy horizontal lines, as in Fig. 
[1] represent Viowfc, and the dashed horizontal fines represent 
the 0/^/3 /3-decay operator in Eq. |T]). 
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FIG. 3. Diagrams in the expansion of the effective interaction defining the two-body part of the second- and third-order Q-box. 
The wavy lines represent Viowfe- We obtain the first- and second-order X-box — the set of all unfolded first- and second-order 
diagrams for the two-body effective operator (not including norm diagrams) — by replacing one interaction in each of these 
diagrams by a /3/3-decay operator (in all possible ways) and restricting the sums over nucleons in the intermediate states to 
either neutrons or protons, as in the first-order X-box diagrams in Fig. |2| 



cates keeping the denominator and numerator separate, 
at the price of introducing disconnected diagrams that 
only cancel when the sum is carried out completely. Here, 
though we evaluate the Q-box to third order and the X- 
box to second order in the interaction, we include only 
one fold in each of the three factors on the left hand side 
of Eq. ([6]), and so opt to follow Refs. [21 [57] in expanding 
the denominator and folding with the numerator. The re- 
sulting expression for the matrix elements of an operator 
A^off is approximatcljQ 

{cd\ Mcfi \ab) = 



where e is the unperturbed energy of both the initial and 
final states (we take the energies to be the same). Both 
Q and X are matrices, with indices corresponding to the 
possible two-body states in the valence space (e.g., a,b 
or c, d in Figs. [I] and [2|. In this paper we report results 
of just the terms explicitly given above, which contain 
between zero and five folds (there is a fold at every matrix 
multiplication) . The terms indicated by ellipses are more 
complicated and presumably less important; they await 
future investigation. 
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^ Because off the need for a special basis, this expression is only 
strictly correct when the terms in square brackets are diagonal. 
They are close to diagonal in the calculations presented here. 



D. Evaluation of Q- and X-Box Diagrams 

We turn now to the Q- and X-boxes themselves, con- 
structed from unfolded diagrams, that we use in Eq. ([7|. 
To construct the Q-box, we take all unfolded diagrams to 
third order in Viowfc- The diagrams appear in Appendix 
A. 2 of Ref. 16 , and the two-body pieces are reproduced 
in Fig. [sj Our X-box has too many diagrams to dis- 
play here, so we characterize the set as follows: we take 
all two-body Q-box diagrams in Fig. p[ and replace one 
interaction line in each diagram (in all possible ways) 
by a /?/3-decay line. We then determine whether each 
intermediate-state nucleon line should be a proton or 
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FIG. 4. (Color online) A Pauli-forbidden two-body diagram 
with a spectator neutron and a three-body diagram, obtained 
by exchanging two ingoing neutron lines, that cancels it ex- 
actly. 



neutron. The result is three times as many X-box di- 
agrams (at second order in Viowfc) as Q-box diagrams in 
Fig.i 

We make one nonstandard choice in evaluating the X- 
box: we restrict the particle lines in the intermediate 
states to be essentially unoccupied. For example, in the 
/3/3 decay of ^^Ge, we omit all contributions from inter- 
mediate protons in the IP3/2 orbit and neutrons in 1^3/2, 
lPi/2: or 0/5/2 orbits, ad we multiply the contributions 
of graphs with intermediate neutrons in the 0179/2 orbit 
by 0.4, its average occupancy. In the decay of ®^Se, we 
omit the same contributions as in ^^Ge and multiply the 
contributions of graphs with intermediate neutrons in the 
0^9/2 orbit by 0.2 and those with intermediate protons 
in the O/5/2 orbit by 0.5. The reason for all this is that 
in a nucleus with more than two valence nucleons, the 
diagram on the left of Fig. [4] — a two-body contribu- 
tion to the /3/3-decay operator with a spectator neutron 
— would be canceled by the three-body diagram on the 
right if we were to include it. By omitting the two-body 
diagrams with intermediate particles in occupied orbits 
we are effectively adding particular three-body diagrams 
(like those on the right of Fig. El) to our calculation. We 
are not including all three-body diagrams, just those that 
cancel Pauli-forbidden two-body diagrams. 

We call this approach nonstandard because it is not 
usually followed in derivations of effective interactions. 
The reason is that in excluding some Pauli-forbidden di- 
agrams, one effectively includes unlinked one- and two- 
body diagrams (see, e.g., Fig. 10 of Ref. [33]) as well 
as the exclusion-enforcing three-body diagrams we want. 
This problem, however, is more pronounced in the Q-box 
than the X-hoK since the latter has no one-body part and 
far fewer ways to unlink diagrams by exchanging lines 
(the horizontal /3/3-decay lines are restricted to have in- 
coming neutrons and outgoing protons). We therefore 
effectively include only very few unlinked diagrams by 
introducing our restrictions in the X-box; the compen- 
sating benefit is a much better account of Pauli exclusion, 
an important physical effect. Diagrams such as the one 
on the left of Fig. |4] result in large contributions that 
should not be present in a full calculation. We cancel 
them with the implicit assumption that the canceling 
contribution from the figure on the right-hand side is 
significantly greater than that of typical third-order dia- 
grams, which we omit. Eventually, though, this assump- 





^''Ge 


«2Se 


Bare matrix element Moi, 


3.12 


2.73 


First-order X-box, without 3p-lh 


5.44 


4.86 


Full first-order X-box 


2.20 


2.40 


First order folded 


3.11 


2.79 


Full second-order X-box 


4.14 


3.92 


Final matrix element 


3.77 


3.62 


CD-Bonn G-matrix 


3.62 


3.45 


N^LO G-matrix 


3.48 


3.33 



TABLE I. The Oiy/Sfi-decay matrix elements Mo„ for ^'^Ge and 
**'^Se at various approximations in our many-body framework. 



tion will have to be tested explicitly. 



III. RESULTS 

To obtain our final corrected shell-model Oi^/3/3-decay 
matrix elements, we combine the individual two-body 
matrix elements of our effective operator with two-body 
shell-model transition densities. Since our aim is a consis- 
tent calculation without empirical adjustment, we really 
ought to take two-body densities from the diagonaliza- 
tion of a valence-space interaction that is derived directly 
from NN+3N forces. While work in this direction is in 
progress, the computation is not yet possible in nuclei 
this heavy. Instead we use two-body densities from ex- 
isting shell-model calculations, the interactions for which 
have been tweaked to fit experimental data in nearby nu- 
clei. For ^^Ge we use the calculation of Horoi [30] and 
for ^^Se that of Ref. [2]; the authors of both have kindly 
supplied us with their transition densities. 

Table H] presents our matrix elements at various levels 
of X-hoK and folding approximations, using Viow k and 
taking intermediate-state excitations to IShuj. Despite 
differences in the NN interaction and size of the basis 
space, contributions from first-order diagrams in both 
^^Ge and ^^Se largely agree with those first identified in 
Ref. [3: particle-particle and hole-hole ladders together 
enhance the matrix element, while the three-particle one- 
hole diagrams cause a dramatic reduction. When folding 
is included, however, the net correction from first-order 
Q- and X-boxes essentially disappears. Taking the com- 
plete set of second-order diagrams into account, we find a 
significant enhancement followed by a modest quenching 
from folding. The final matrix element is approximately 
20% percent larger than the bare matrix element in ^^Ge 
and about 30% larger in *^Se. The primary reason for 
the different effects in ^^Se and ^^Ge is the difference in 
the omitted intermediate-state orbits discussed in Sec- 



tion IIP If we include those orbits, as is standard prac- 
tice in the construction of effective interactions, the ma- 
trix element is reduced by about 10% in ^^Ge and 15% in 
^^Se. In Ref. [llj . which contains a preliminary account 
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TABLE II. Convergence of 0!//9/3-decay matrix element in *^Se 
with respect to allowed intermediate-state excitations. In all 
cases we work in a harmonic-oscillator basis of 13 major shells. 



of our calculations in *^Se, we obtained 3.56 instead of 
3.62. The small difference is due to the inclusion in Ref. 
[llj of Q-box restrictions and the addition here of a term 
in the expansion of the norm denominator. Though the 
two results are close, we believe that the one reported 
here is likely closer to the real matrix element. 

Several other aspects of the calculation are robust. As 
seen in Table |llj our Viowfe results at ISHu are well con- 
verged to 3 or 4 digits. And as Table [l] shows, chang- 
ing the interaction to a G-matrix (in 11 major oscillator 
shells) in place of Mow/c does not affect the results sub- 
stantially. Finally, although we emphasized our proce- 
dure of requiring intermediate-particle lines in X-box di- 
agrams to be unoccupied in the nucleus in question, other 
prescriptions yield similar results once norms and folding 
are included: in *^Se, for example, we obtain a final ma- 
trix element of 3.50 if we restrict particle lines in both 
the Q and X boxes, and 3.03 if we impose no restrictions 
at all. We should note, however, that at various interme- 
diate stages of the calculation, the procedures yield quite 
different results. And other parts of the calculation leave 
room for change as more physics is included. 

We turn now to a discussion of 2^/3/3 decay. As noted 
above, we use the closure matrix element M|J, as a proxy 
for the full matri x ele ment, a step that limits how much 
we can say. Table III shows the matrix element for ^^Ge 



with and without the intermediate-state restrictions we 
impose on occupied or partially occupied orbits in calcu- 
lating Mqi,. Imposing the restrictions here increases the 
matrix element, as in QvPfi decay; in this case, however, 
the result is probably undesirable, given that shell model 
calculations of 2z^/3/3 decay in ^^Ge overestimate M2y- 
On the other hand, omitting the restrictions increases 



Restricted Unrestricted 

Bare matrix element M2J, 0.57 0.57 

First-order X-box, without 3p-lh 0.99 0.89 

FuU first-order X-box 0.37 -0.60 

Full second-order X-box 0.79 -0.37 



Final matrix element 



0.96 



0.70 



TABLE III. The 2!//3/3-decay closure matrix M|i for '^'^Ge at 
several levels of approximation, with and without restrictions 
on occupied intermediate-particle lines. 





FIG. 5. (Color online) Schematic representation of diagrams 
contributing to renormalization of qa in 2i//3/3 decay. 



the negative contribution of the 3p-lh diagram to such 
an extent that the matrix element changes sign. The sign 
is eventually reversed by higher-order contributions and 
folding, ultimately yielding a result that is approximately 
unchanged from the bare matrix element. It is difficult 
to be comfortable, however, with a low-order correction 
that changes the sign of the matrix element. The sensi- 
tivity of the numbers suggest that terms with more folds, 
of higher order, or involving more valence orbitals could 
also have a significant effect. 

Another reason (aside from the sign changes in the 
right-hand column of Table III I for preferring to restrict 



intermediate-sate orbits in the X-box is connected to the 
long-standing problem of the apparent suppression of the 
axial-vector coupling constant qa in the nuclear medium 
[21]. While the suppression probably has many sources, 
configurations outside the valence space are likely to play 
a key role. Though the bare operator governing weak de- 
cay is one-body, we can simulate the effect of qa suppres- 
sion in 2uj5j3 decay by including only closure diagrams 
that have the form shown in Fig. [5] Such diagrams, in 
which only a single 2z^/3/3-decay line connects the two 
nucleons, incorporates only the renormalization of the 
one-body weak current. 

When we base our calculation of M|J, on only these 
diagrams and at the same time account for occupied 
intermediate-state orbits, we find in *^Se that the full 
result is smaller than the bare result by 38%, implying 
an effective qa of about 1.0, a reasonable value (in ^^Ge 
the value is about 0.7). On the other hand, when we 
take no account of occupied orbits we find (again, with 
only diagrams of the form in Fig. [s] included) that the 
closure matrix element changes sign, something that is 
impossible through the renormalization of qa alone. The 
sign change reflects the same strong effect of the 3p-lh 
diagrams observed in Table [Till Though these consider- 
ations are not conclusive, they do indicate that taking 
the Pauli principle into account is a beneficial. Work is 
currently underway to investigate qa quenching more di- 
rectly. By focusing on the one-body operator, we can use 
analogous many-body techniques, based again on chiral 
NN and 3N physics, with the effects of two-body cur- 
rents implemented consistently in the bare operator |32| 
to understand the origin of gA quenching. 

Whatever the outcome of that investigation, it is clear 



that the 2i//3/3 matrix element is sensitive to many details 
in the wavefunctions, much more sensitive than its 0i>/3/3 
counterpart. Thus, although the increase of the 2i//3/3 
closure matrix element does not bolster the case for our 
QvjSp calculation, neither, in our view, does it weaken it 
much. 



IV. DISCUSSION AND OUTLOOK 

We have used chiral nuclear forces and many-body 
perturbation theory to calculate an effective shell-model 
Qv j3 j3-dec&y operator, taking into account corrections to 
the bare operator from configurations outside the valence 
to second order in the interaction. The resulting nu- 
clear matrix element is approximately 20% larger than 
the bare matrix element in ^^Ge and about 30% larger 
in ®^Se. These new results represent our current best 
estimates for the matrix elements but probably do not 
tell the whole story. We have omitted a number of ef- 
fects that could further alter the results. To do better, 
we must first establish consistency between the Hamil- 
tonian and our effective operator. This will require the 
construction of full non-empirical valence-space interac- 
tions in the pf shell from NN and 3N forces; work in that 
direction is in progress. A related improvement will be 
to include the effects of chiral 3N forces in the X-box, in 
addition to the effects of chiral two-body currents in the 
bare operator [32- 

At the many-body level, the importance of third- and 
higher-order terms in the X-box and additional folding 
contributions must be understood. Since we have found 
the effects of bubble diagrams to be the most important 
in our perturbative expansion, it would be worthwhile to 



pursue a nonperturbative calculation of the effects of core 
polarization (which these diagrams represent), like that 
done for effective interactions in Ref. [S]. Perhaps the 
most significant obstacle to a truly reliable result, how- 
ever, is the implementation of induced three-body oper- 
ators. Recent work j33l[34] indicates that such operators 
are not negligible, and even here we have shown that 
three-body diagrams of the form in Fig. [4]are important. 
Unfortunately, the number of induced three-body dia- 
grams is so large that nobody has computed them even 
in the construction of effective interactions. We must 
find a way to at least estimate their size if we want to 
pursue perturbation theory to its conclusion. Controlled 
nonperturbative approaches [7l[9j are on the horizon, but 
the inclusion of induced three-body terms is technically 
difficult there as well. In none of these approaches is the 
problem impossible to overcome, but doing so will require 
diligence and creativity. 
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